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SUMMARY 

A total strainrange version of stralnrange partitioning (SRP) was 
proposed by Halford and Saltsraan to enhance the manner in which SRP is applied 
to life prediction. This report describes, for the SRP model, how advanced 
reliability technology can be applied to (a) perform risk analysis and (b) 
to derive safety check expressions. 

Uncertainties existing in the design factors associated with life 
prediction of a component which experiences the combined effects of creep 
and fatigue can be identified; (a) inherent uncertainty in material behavior, 
(b) statistical uncertainty associated wi;h parameter estimates resulting 
from small samples of fatigue specimens, (c) modelling error associated 
with the SRP model, (d) data scatter in the environment, e.g., loads, 
temperatures, hold times, (e) modelling error associated with service strain 

analysis. Examples are presented which illustrate how reliability 
analyses of such a component can be performed when all design factors in 
the SRP model are random variables reflecting these uncertainties. 

Using the Rackwitz-Fiessler and Wu algorithms, estimates of the safety 
index S and the probability of failure p^ are demonstrated for an 
SRP problem. Methods of analysis of creep-fatigue data with emphasis on 
procedures for producing synoptic statistics are presented. An attempt 
was made to demonstrate the importance of the contribution of the uncertainties 
associated with small sample sizes (fatigue data) to risk estimates. In 
the example presented, the influence of such statistical uncertainty was 


small. 
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Finally, an Illustration of the procedure for deriving a safety check 
expression for possible use in a design criteria document was presented. 
The format employs partial safety factors (PSF) which are derived from re- 
liability analyses. The safety check inequality has the appearance of a 
"conventional” design requirement, and therefore is familiar to designers. 
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NOTATION 

Constant in linear model; defined by Eq. A. 2 
Least squares estimator of a 
Coefficients; Eq, D.6 

Coefficient of strain-life curve; See Eq. A,1 
Coefficient of strain-life curves; defined in Eq, 6 
Median of A^ 

Constant in linear model; defined by Eq, A. 2; also exponent in 
inelastic strain-'life curve 

Least squares estimator of b 

Coefficient of elastic strain-life curve in which only PP strain 
is present; defined by Eq. 10 

Median value of Bp^ 

1/b 

Hysteresis loop in which tensile creep reversed by compressive creep 
Hysteresis loop in which tensile creep reversed by compressive plasticity 
COV of A^ 

Coefficient of variation 

Exponent of elastic strain-life curves; defined by Eqs. 10 and 11 

Coefficient of B - B relationship; defined in Eq. 12 
PF i 

Type I extreme value distribution of maxima 

Fraction of the total of each strain range type i = PP, PC, CP, CC 
Probability density function of 

Probability density function of ACg 
Mean value of fpp 
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NOTATION - (continued) 


fp^ Mean value of fp^ 

g Function which accounts for statistical scatter; defined in 

Eq. A. 8; also Eq. £.3 

G A random variable which quantifies modelling error in computing 

service strain range 

H A random variable which quantifies material behavior uncertainties 

in computing service strain range 

J Coefficient of inelastic strain-life relationship; defined by Eqs. 8 & 9 

LN Lognormal distribution 

n Sample size 

N Cycles to failure; also normal distribution 

N Service life 

o 

N^ Cycles to failure for i strain range type; i = PP, CP, PC, CC 

Probability of failure in service life 
Target risk or probability of failure 
P(‘) Probability of 

PC Hysteresis loop in which tensile plasticity reversed by compressive 

creep 

PP Hysteresis loop in which tensile plasticity reversed by compressive 

plasticity 

PSF Partial safety factors 

Q Load (or nominal stress) range on the component 

s Sample standard deviation; estimate of a 

SRP Strain range partitioning 

t Time 

t , Students* t variate 
a;n-l 

T Hold time; in general a random variable 

T Mean value of T 
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NOTATION (continued) 

Reduced, coordinate; defined in Eq. D,3 
X Los^qAe 

til 

X^ Log^pAe^ where i refers to i specimen 

Y a + bx; defined by Eq. A. 2; also log^^N 

A 

Y Least squares line 

th 

^i lO^i ^ refers to i specimen 

z Standard normal variate 

a 

* As a superscript, refers to design point for that variable 

a Reference level for g; See Eqs. A. 7 and A. 8 and Refs, 5 and 11 

3 Safety index 

3^ Target safety index 

Y Empirical function of 0 and T; defined by Eq. 2 and Fig. 5 

y Partial safety factor for variable X 

A 

6 Exponent of Bpp ~ relationship; defined in Eq. 12 

Ae Total strain range 

Ae^ Strain range; i = PP, GP, PC, CC 

inelastic strain range 

ACpp PP strain range 

Ae Cr strain range 

Cr 

ASpc PC strain range 

Ae^^ CC strain range 

ACg Total service strain range 

Ae„ Mean value of Ae_ 

b b 


NOTATION - (continued) 


Total strain range to produce failure at life N; describes the 
strength of the material; defined in Eq. 14 

Temperature; in general a random variable 

Empirical function of 0 and T; defined by Eq. 2 and Fig. 5 

Equivalent normal mean 

Standard deviation of y|x 

Equivalent normal standard deviation 

Equivalent standard deviation 

Standard normal density function; also empirical constant 
defined by Eq. 2 and Fig. 5 

Standard normal cumulative distribution function 



TABLE OF COMPARISONS OF NOTATION 


The authors of this paper used notation which differs from that of the 
authors (Halford and Saltsman [ 2 ]) of the advanced version of SRP. The 
intent of this change of notation was to simplify the presentation in a 
reliability format. Following is a comparison of notation of some key 
parameters 


Symbol 


Halford 

Saltsman 

This 

Paper 


t 

T 

hold time 

Ae , 
el 

Ae 

e 

elastic strain range 

^^in 

Ae , 
in 

inelastic strain range 


a^t 

total Strain range 


J 

coefficient of the inelastic strain life curve 

c 

c 

exponent of the inelastic strain life curve 


B. 

1 ' 

1 

j 

coefficient of the elastic strain life curve for the ; 
strain type; i = PP, PC, CP, CC 

b 1 

d 1 

exponent of the elastic strain life curve 


f . 

1 

th 

fraction of i strain type to total inelastic strain 


D. 

1 

coefficient of relationship between Bpp and 

a 

6 

exponent to hold time in relationship between B„ 
and B_ 


1. INTRODUCTORY COMMENTS ON STRAIN-RANGE PARTITIONING 

The method of strainrange partitioning (SRP) for predicting high 
temperature low cycle fatigue was introduced a decade ago by Manson, Halford, 
and Hirscliberg [1], This scheme for making life predictions is based on 
explicit knowledge of the magnitudes of ],he inelastic creep and plastic 
strains present in a cycle of loading. Unfortunately for typical engineering 
applications, the magnitudes of the plastic strains are small and they 
cannot be calculated reliably fror.! nonlinear structural analysis methods. 

Halford and Saltsman have proposed a method which enhances the manner 
in which SRP is applied to life prediction [2}. They developed the basic 
Manson-Coffin plastic strain-range power law of low cycle fatigue into a 
total strain-range representation by the addition of the elastic and plastic 
strain life relationships. It is argued that this method, a total strain- 
range version of SRP, has the promise of more accurately estimating cyclic 
lifetimes over a much broader range of strains and lives than was possible 
on the basis of either the plastic or elastic strain-range alone. 

Many uncertainties exist in the process of employing SRP for life 
prediction. In a broad sense, these would include (a) scatter in environ- 
mental data, and uncertainty in the computations of the environment, e.g., 
temperature, (b) modelling error associated with the procedures for comput- 
ing loads on the components and then computing responses (stresses), (c) 
uncertainty in the responses of the material to the environment, (d) scatter 
in fatigue data, (e) modelling error of the theoretical strength model, 
i.e., SRP. The general goal of this study is to demonstrate how modern 
probabilistic design theory can be employed to predict reliabilities of com- 
ponents subjected to high temperature low cycle fatigue. SRP will be the 
basic prediction method used. 
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For reference purposes, the following basic definitions and descriptions 
of SRP are Included from Ref. 3 . First consider a hysteresis loop as 
shown in Fig. 1. Defined are the inelastic elastic and total 

(At) strainranges . The basic premise for SRP is that in any hysteresis 
loop there are combinations of just two directions of f»training and two types 
of ivielastic strain. The two directions are, of course, tension (associated 
with a positive inelastic strain rate) and compression (associated with a 
negative inelastic strain rate); the two types of inelastic strain are time 
dependent (creep) and time independent (plastic) . It should be noted that 
only a portion of transient creep strain should be considered as plastic strain 
and only the steady-state component be considered as creep strain [A]. 3y com- 
bining the two directions with the two types of strain, we arrive at four possible 
kinds of strainranges that may be used as basic building blocks for any conceivable 
hysteresis loop. These define the manner in which a tensile component cf 
strain is balanced by a compressive component to close a hysteresis loop. 

The types of strain are illustrated in Fig. 2 and are described as follows; 

(a) Tensile plasticity reversed by compressive plasticity is designated 
a PP stralnrange and represented by ASpp. 

(b) Tensile creep reversed by compressive plasticity is designated a 
CP stralnrange and represented by Ae__, 

ur 

(c) Tensile plasticity reversed by compressive creep is designated a 
PC stralnrange and represented by ACp^. 

(d) Tensile creep reversed by compressive creep is designated a CC 
stralnrange and represented by A6__. 


CREP 


PLASTIC aow 


INELASTIC STRAINRANGE.Ag.j^ ELASTIC STRAINRANGE, Ac 

TOTAL STRAINRANGE, Ac 


Fig. 1. Hystertris Loop, 
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Fig . 2 . Idealized hysteresis loops for the four basic types of 
inelastic strainrange. 
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The notation lor the subscripcs for the stralnranges uses the type of 
tensile strain first, fe.Uowed by the type of compressive strain. The name 
strainrange partitioning was chosen because it represented the premise that, 
in order to handle a complex high-temperature, low-cycle fatigue problem, 
the Inelastic strainrange must first be partitioned into its components. 

The strength of the material is described by strain life curves, 
an example of which is shown in Fig. 3. These relationships follow the basic 
Manson-Coffin law. 

Given a stable hysteresis loop under constant amplitude oscillatory 
loading, as shown in Fig. 1, the fraction of each strainrange typ? f ^ , a 

td, 

component of the total Inelastic strainrange, is identified using an 

algorithm as des irlbed in Refs. 1 and 3 , For example, f_„ = AG__,/Ae. . 

PP PP in 

4 4 

' I Ae^ I ^1 “ 1 i - PP, CP, PC, CC (1) 

i=l i~l 

Finally, it should be noted that notation of the original SRP work has 
been changed somewhat herein. This was done for mathematical convenience 
in applying reliability theory, 

2. UNCERTAINTIES IN THE LIFE PREDICTION PROCESS 

For typical designs in a high temperature environment, the present 
state-of-the-art precludes an accurate deterministic definition of the €>n- 
vironments and associated material responses. Moreover, fatigue behavior 
under carefully controlled conditions is characterized by significant 
uncertainty as evidenced by the large scatter in fatigue failure data. 

The goal of this study is to c?..st the total strain range version of 
SRP into a reliability format. All sources of uncertainty will be identified. 
Techniques for quantifying uncertainty will be addressed. Mechanisms for 
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formal introduction into the limit state function will be described as will 
modern methods for performing the reliability analyses. 

The sources of uncertainty in the process of fatigue life estimation 
can be identified as follows; 

1. Environment. There is uncertainty in temperatures, hold times 

and pressures, static and centrifugal loads, etc. Stress producing environ~ 
ments may be random processes or deterministic processes with random magni- 
tudes, Statistical descriptions of the environment may be available. 

2. Response to the Environment. The computational methods for computing 
stresses will contain modelling error resulting from the assumptions made. 

3. In general, strain at the fatigue critical point will be a random 
process. Strain range and mean strain for each hysteresis loop will be a 
random variable, reflecting uncertainties in material properties (e.g. , 

Young’s modulus) as well as environmental and analysis uncertainties. 

4. Dividing hysteresis loops into strain types. The process of 
Identifying the fractions of each hysteresis loop associated with each strain 
type (PP, PC, CP, CC) will likely contain uncertainty. The method used may 
not accurately reflect real strain behavior of the material. 

5. Linear damage rule. The interaction damage rule used in 

SRP may not accurately describe fatigue behavior of the material [ 5 ] • 

This uncertainty is referred to as modelling error, and is associated with 
the theoretical model which is assumed to define strength. 

6. Material behavior. Fatigue data is typically characterized by 
’’large" scatter. Moreover, parameters used to provide statistical summaries 
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of the data are themselves random variables when the estimators are used 
to represent the parameters. 

7. Material behavior, . . . other uncertainties. The fatigue strength 
of a material may be influenced by processing operations, e.g., cold work- 
ing and heat treating), . . . and assembly operations (e.g,, bolting, 
shrink fits). Uncertainties in material strength may result. Moreover, 
material strength may be influenced by time and/or by corrosion and/or 
extreme thermal environments to a degree which is not accurately known. 

Following are discussions of the components of the SRP model and a 
demonstration of how modern methods can be used to perform reliability analysis 
on a high temperature low cycle fatigue problem. 

3. SERVICE STRAIN 

It is assumed that the component operates at constant temperature and 
that the temperature is high enough so that creep deformation must be con- 
sidered. Also, it is assumed that the loading is constant amplitude. The 
physical problem is illustrated in Fig, 4. Assume that the load (or nominal 
stress) range, Q, is a random variable reflecting (a) uncertainties in the 
environment, and (b) modelling error in translating the effect of environ- 
ment to loads on the component. 

The total service strainrange Ae at the notch will be a function of 

D 

Q, temperature 0 and hold time T as shown in Fig. 5. An analytical model for 
service strain range can be formulated as 

ASg = yQ + (<I>Q)’^ 


( 2 ) 
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Fig. 5. Relationship Between Service Load Ranges and Notch Strain 
Range. 
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where the parameters y» and p may be functiuns of o and T. 

Y “ y(0j1) 

<!> » <J>(0,T) (3) 

n » n(C,T) 

Both hold time and temperature can be random variables by virtue of uncertainties 
in the operating environment and perhaps the codes used for their prediction. 

Thus Y» '!>) and n will, in general, be r?ndom variables. 

There are two other sources of uncertainty here. First, the method by 
which strains are computed from load will contain modelling error. Then 
there will be uncertainties in the material response as could be measured 
from experimental data. A more general form of Ae would be 

AEg = GH (yQ + (<})Q)'^) (4) 

where G and H are random variables which account for modelling error and 
material behavior respectively. 

4. IDENTIFYING THE STRAINRANGE COMPONENTS 

The SRP literature describes the mechanical procedure for quantifying 
the partitioned strainrange components of a complex hysteresis loop [ 3 ] . 

It is possible to have only three of the four types in the same loop. 

Let fjj^; i = PP, PC, CP, CC denote the fractions of each partitioned 
strainrange. The sum of the fractions is unity and as an example 
consider 
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^ ' *PP ^PC ^ *CC 

Each term can be considered as a random variable. Three sources of uncertainty 
can be identified, First, there may be uncertainty in the way that the loop 
is analyzed (modelling error) or there could be some error in the basic 
algorithm for dividing the plastic strains. It is expected that this error 
may be small and difficult to quantify. Nevertheless, if the can be 

modelled as random variables, no problem is presented to the reliability 
method. 

But the fj^’s will also depend upon hold time T and temperature G. 

If T and 0 are known random variables, and if their functional relationship 
to f^ can be described, then in theory, the distribution of each f^ can be 
derived. Fig. 6 illustrates the relationship which must be established from 
testing. 

In the design equation, the f 's are clearly not independent as seen 

i 

from Eq, 2, For three strainrange types, two f^'s can be specified indepen- 
dently, and the third f^ expressed as a function of the other two. A demon- 
stration of how to handle this in a reliability format is provided in the 
examples below. 

5. THE STRAIN-LIFE RELATIONSHIPS. HOW SCATTER IN FATIGUE DATA IS TREATED 
The inelastic strainrahge-life curves are established by conventional 
SRP techniques. It is assumed that the data will follow a linear trend on 
log- log paper (the Manson-Cof fin law) and that the techniques of basic linear 
model analysis apply. Methods of analysis for e-N data for design purposes 
are summarized in Appendix A. 
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Ab shovm in Fig. 7, it is also assumed that the slope, b, of each 
e~N curve is the same and equal to the PP curve (for which the sample sizes 
are usually much larger). Fig. 3 illustrates how this assumption may be in 
error, but these curves were based on a small mimber of points. 

The empirical relationship for each strainrange type is given as 


N * A (Ac 


^PG “ 


^CP “ 


^cc “ 


( 6 ) 


The e-N curves of Fig, 7 are median curves through the data and are 
defined by the relationships given on the figure. The tildes indicate 
median value. 

Scatter in observed fatigue data is accounted for by treating the A's 
as random lognormally distributed variables. The exponent b is considered 
to be constant. Appendices A and B describe the process of translating 
e-N data into statistical parameters of the random variables A^. 

In order to construct the appropriate inelastic strain-life curve, the 
basic SRP model is employed [1]. The total cycles to failure, N, is 


4 

i f /N i = PP, PC, CP, CC (7) 

i=l ^ ^ 


4 



Rearranging, the resulting e-N relationship becomes. 


Inelastic 

Strain 

Range 
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whe'i^e, 




J 



1 c 


i 


c “ 1/b 


( 8 ) 


(9) 


Eq. 8 then provides an expression of the inelastic strain life curve, An 
example is given in Fig, 8 in which it is assumed that only PP and PC strain 
are present, 

The elastic strain range-life curves are ertablished as follows. 

First, the PP line is defined from the data using the method of Appendix A. 
The strain life model is, 

Ae = (10) 

e rr 

in which B^p Is a random variable and the exponent d is a constant. This 
curve is illustrated in Fig. 9. Then the elastic strain life curve for a 
given hold time and for a given constant amplitude load can be established 
from experimental data in the same way. 

As the hold time is increased, the e-N curves will Indicate lower 
fatigue strength as suggested by Fig. 9. The curves will be parallel to the 
PP curve and will have the form 

Ae -= B.N^ i « PC, CP, CC or (11) 

e i 

some combination 
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F1g. 9, The Elastic Strain Life Curves. 

(Note: It is possible for another line to be above the 
PP line, e.g., for a material which experiences strain 
age hardening.) 
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The subscript i refers to the form of the hysteresis loop, i.c,, PC, CP, CC, 
or a combination of these with PP, 

Experimental data has suggested that the assumption of (a) parallel 
e-N lines and (b) an empirical form 

- D.t’’ (12) 

are reasonable. An illustration of data which supports Eq. 12 is shown 
later in Fig. 16, In general, and 6 will be functions of hold time, 
temperature, and type of cycle, Using the general scheme as de- 
scribed in Appendix A, will be a random variable, and 6 will be constant, 
and 6 will be established from experimental data so that the distribution 
of will reflect both material and statistical uncertainty. 

In general, the elastic strain-life expression is given by Eq. 11, 
substituting from Eq. 12, 

= Bpp expl-Djl'^l (13) 

Combining the plastic and elastic strain ranges, the total strain 
life curve is given as 

= Ae + + JN^ (lA) 

i e in i 


The strain-life curves are illustrated in Fig. 10. 

Upon substituting the expressions for J and B^. the total strain life 
expression becomes 


Ae. 


{B 


PP 


exp [ “D^T 


1 ) + 


{Jf.M }V 
j ^ ^ 


(15) 
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This is the definition of fatigue strength. It is assumed that B^p, 
and are random variables. Therefore, for a given life 
also will be a random variable. 

6 . RELIABILITY A1>!.M.YSIS 

The probability density function (pdf) of the fatigue strength, denoted 

as f , will be a function of N; it is illustrated in Fig. 11 at the intended set 

Etji 

vice life N , The service strain range is denoted as Ae„. Also shown on this 

O a 

figure is the pdf of Ae„, denoted as f . It is assumed that the strain 

S Cg 

range will be constant over the life of the component, but the magnitude 
(ACg) is treated as a randoBi variable to reflect uncertainties in the 
environment as well as the procedures used to compute the strains. 

The event of failure is defined as (As„ < Ae„) , and the probability 

X D 

of failure is 

p^ = P(Ae,j^ < ACg) (16) 

In the language of mechanical reliability, Ae.j, is the "strength,” and 
ASg is "stress." 

EXAMPLE 1 

Consider a component, subjected to a constant amplitude oscillatory 
stress, which is expected to experience some inelastic strain, of the PP 
and PC types only. Thus, the fatigue strength of the material would be de- 
fined by a special case of Eq. 15, 

f f 

AEt = exp[-Dp^ T*^]} n'^ + (17) 

which includes only PP and PC strain terms. 
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Fig. 11. Relationship Between the Distribution of 
Service Strain and Fatigue Strength of 
the Material 
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The material is to be a nickel base alloy AF2-1DA at 1400 F (760 C) , A 
summary of the mechanical properties is given in Table 1. This example il- 
lustrates the use of modern reliability methods to assess the structural 
performance of the component. For a specified life, and the associated 
service strainrange, it is required to estimate the safety index and the 
probability of failure. Data used for all design factors are summarized 
in Table 2. Commentary on how the parameters are determined from these 
data is provided in the following. 

Hold tim e, fpp, fp^. Assume that there is uncertainty in the hold time, T. 
Thus, T can be a random variable, and for this example the coefficient of variatic 
was assumed to be only 5%, But fpp and fp^^ will be functions of T, and it is 
therefore necessary to provide explicit functions. Fig^. 12 shows how such a re- 
lationship might appear. This is used for the example and is not based on actual 
data. In fact, such a relationship could be established from a simple test. 
Scatter in material behavior is not considered here. It should be noted 
that the f__~f_ -T relationship will also depend upon temperature 0 and 

rr r\j 

service strain range, Ae ; the latter is also a function of T and G. A 
simplified physical model is employed herein for demonstration purposes. 

In this example, it is assumed that the uncertainty in hold time 
will be relatively small, and that in the first approximation the relation- 
ship between fpp and T is the tangent to the curve at the mean value (in 
this example, 100 sec.) as shown in Fig. 12. 

fpj, = 1.1 - 0.10 log^Q T (18) 



( 19 ) 


Thus, upon substituting Eqs. 18 and 19 into Eq. 17, the plastic strain 
life curve can be expressed as a function of T. 
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TABLE 1 

Mechanical Properties of AF2-1DA 

Young’s Modulus 
Yield Strength 
Ultimate Strength 
Reduction of Area 

Stress Rupture Properties 

Stress (ksi) Reduction 
in Area (%) 

135 15.8 

130 14.6 

125 15.0 


at 1400 F (760C); Ref. 6 

25 X 10^ ksi 
123 ksi 

164 ksi i 

22.3% 

Time to I 

Rupture (hrs) ; 

1.1 i 

. I 

2.1 j 

196. '! 

f| 


1 

1 

' I 


! 

i . I 

I i 
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TABLE 2 

Data for Example 1 


Random Variables 

Distribution 

Mean 

> 

o 

<_> 

A&S 

EVD 

4.45E-3 

20 

T 

N 

100 

5 

Bpp 

LN 

0.0216(b) 

9.9 

°PC 

LN 

0.0447(b) 

30 

fpp 

N 

0.80 

5 


LN 

0.0281(b) 

47 

A2 

LN 

0.0156(b) 

69 


Constants 


1 

' 0.25 

II 

O 

-0.637 

d 

-0.117 

Nq, service life 

1 ,000 cycles 


Notes 

(a) Abbreviations 

EVD Type I extreme value distribution of maxima 
N Normal 

LN Lognormal 

(b) For lognormal variates, the median is used rather than the mean 
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Data used for all of the design factors are summarized In Table 1. 
Commentary on these data Is provided In the following, 

Stress . It Is assumed that the service strain amplitude, ASg, has 
a Type I extreme value distribution of maxima, (EVD) . The COV of 20% is fairly 
typical of loading variables and reflects primarily modelling error resulting 
from assumptions made in the computational procedures which translate en- 
vironment into notch strains. A more refined and complete model in which 
ASg is expressed as a function of temperature and hold time is presented 
above in Eq. 4, but the simplified approach is used here simply to illus- 
trate the reliability methods. 

FP Strain range Data , The PP strain range-life data for AF2-1DA at 
1400 F is shown in Fig. 13. Methods of analysis of these data are described 
in Appendix A, and a summary of the results is given in Table 3. It should 
be noted that the uncertainty in fatigue strength is described by the random 
variable A^, whose COV includes data scatter as well as statistical uncertaintic 
in the estimates of the least squares parameters. 

PC Strain range Data . The PC strain range-life data is shown in Fig. 14. 

To analyze the data, it is first assumed that the slope will be the same as 
the PP curve, i.e., b = - 1.57. A least squares method, with the exponent 
known, is employed, and the results are summarized in Table 4. The COV 
of A^ reflects both data scatter and statistical uncertainty, the latter 
which is quantified using the methods described in Appendix B. 

Elastic Strainrange Data . The PP elastic strain life data is plotted 
in Fig. 15. A summary of the statistical analysis of this data is pro- 
vided in Table 5. Basic analysis methods are summarized in Appendix A. 

Note that for this data, the least squares analysis is applied to the form 
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TABLE 3 

STATISTICAL ANALYSIS OF PP DATA 
(AF2-1DA at I400F, 760C) 

• Transformation 


Yi « log^Q N^ 


“ logj^QAt^ 


• Data (See Fig. 13); Sample Size, n ■ 9 


4Epp(%) 

Nf 

(cycles) 



"i 

.896 

43 


-2.047 

1.633 

.368 

200 


-2.434 

2.301 

.154 

756 


-2.812 

2.878 

,104 

1,322 


-2.983 

3.121 

.089 

2,695 


-3.0506 

3.430 

.037 

4,205 


-3.432 

3.624 

.032 

5,745 


-3.495 

3.759 

.018 

25,433 


-3.745 

4.405 

.011 

59,121 


-3.959 

4.772 


• Least Squares Analysis [See Appendix A] 

A 

Y = a + bx 

a = -1.552 b = -1.57 

s = 0.138 


Statistical Model (See Appendix A for detail of this example) 
^PP ^l^^^PP^ 


b = b = -1.57 


Median of A^; “ 0.0281 

COV of A^; = 0.472 


PC Inelastic Strainrange, ae 


30 


oRiGim p:/ ’ 



Cycles to Failure, N 


Fig. 14. PC Strain Range-Life Data for AF2-1DA at 1400 F. 


3 


TABLE 4 

STATISTICAL ANALYSIS OF PC DATA 


• Transformation 


logio h 




• Data (See Fig. 14); Sample Size, n * 6 


pc(^) 

Nf 

(cycles) 


^1 

.671 

51 


-2.173 

.290 

212 


-2.537 

.184 

300 


-2.735 

.069 

904 


-3.158 

.043 

1,807 


-3.366 

.051 

3,380 


-3.291 


Y 


i 


1.7 

2.3 


2.47 


2.9 

3.2 


3.52 


• Least Squares Analysis (assume the same slope as the PP data; 
b = -1.5?) 

Y = a + bX 


a = -1.807 


s = 0.173 

• Statistical Model (See Appendix B for detail of this example) 

^PC " ^2*''^^PC^ 
b = b = -1.57 
Median of ~ 0.0156 

COV of A.; C, =0.687 
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TABLE 5 


Statistical Analysis of Elastic Strain Lific Data 
(AF2-1UA at 1400F) 


• Transformation 




^1 " 


• Data (See Fig, 15); Sample Size, n = 9 


iEg(X) 

(cycles) 



''i 

1.492 

43 


-1.826 

1.633 

1.155 

200 


-1.937 

2.301 

.898 

756 


-2.047 

2.878 

.898 

1,322 


-2.047 

3.121 

.799 

2,695 


-2.097 

3.430 

.800 

4,205 


-2.097 

3.624 

.783 

5,745 


-2.106 

3.759 

.703 

25,433 


-2.153 

4.405 

.652 

59,121 


-2 . 1B6 

4.772 


• Least Squares Analysis (See Appendix A) 

A » 

Y = a + bx 

a = -14.21 b = -8.532 

s = 0.2581 

• Statistical Model 

N = A(Ae^)^ 
b - -8.532 
Median of A = ^ = 10 

COV of A, = 1.01 (See Eqs. A. 7 ard A. 10; g(.01,9) = 1.41 used 
here) 
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TABLE 5 (continued) 
• Alternate Form Used in Analysis 

i) (1 “ 1/b 


ii) 

PP 

(1/X)^^^ - 0.0216 



[ „ a/b)^ 1^2 

iii) 


^ } (1 + cp -1 = 0.099 


The relationships of ii) and iii) are valid only when A (and therefore, B) 
have lognormal distributions. These are basic forms for lognormal variates 
[e.g. , See Ref. 5] . 
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N “ A(Ae)^. For the strength formulation, statistics on the parameters 
of the equation As = B N ^ are required. Forms for relating statistics 
between A and B anu c and d are given in Table 5. 

Elastic Strain-Life vs. Hold Time Relationship . The relationship 
between elastic strain-life curves and hold time is established from experi- 
mental data as shown in Fig. 16. Bata from CC and CP strains, not shown, 
supported a selection of a slope of 0.25 for the data. Thus, the empirical 
form, relating Bpp and Bp^ is, 

tn(app/Bpg) ■= Dp^T°-“ (20) 

Least squares analysis is performed; the statistics on Dp^^ are presented 
in Table 6. 

Elastic and Inelastic Strainrange-Life Relationships . The strain- 
life curves, employing the total strainrange version of SRP and the data 
of Table 2, are presented in Fig. 17 for reference only. These curves 
suggest that the influence of creep in this example is relatively small. 

Reliability Analysis . The fatigue strength, Ae^, of the material is 
given by Eq. 17 with substitutions of Eqs. 18 and 19 for fpp and fp^,. 

The event of failure is (A^g < Ae^) . The following methods will be used 
to evaluate the probability of failure, 

1) Monte Carlo. This method is widely employed for solving compli- 
cated probability problems. It is a very useful tool, but it suffers 
from high computer costs relative to accuracy. 

2) Rackwitz-Fiessler (R-F) . The R-F scheme is a numerical method for 
evaluating reliabilities in problems such as this one [ 5 ] . It 

is now widely employed and details of the method are well docu- 
mented I 7 , 8 , 9 ]. It has been demonstrated by Wu et al. [10] 



ORIGINAL ® 

OF POOR QUALITY 

Lognn(Bpp/B,.)l 

PC ^ 








CC; Ln(Bpp/Bj,(,) = 0.n84(t)- 
PC: Ln(Bpp/Bp|,) = n.fl46(t)- 

CP: Ln(B„/B.„) = 0.040(t)' 


Log t, (Hold Time, secs.) 


Fig. 16. Time-Dependent Intercepts for Elastic Strain Range Life Relations 
AF2-1DA, 760°C, Halford and Nachtigall Data [6]. 
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TABLE 6 

STATISTICS ON Dp^ 


• Transformation 

X ■= log3,gt 

• Data (See Fig. 16); Sample Size, n = 6 


1 


-1.00 

1.83 

-0.813 ’ 

1.86 

-0.779 1 

2.62 

-0.733 

2.46 

==0.707 : 

2.51 

-0.466 

3.16 


* Least Squares Analysis 
Y = A + 0.25 X 

Least squares relationship 

A = il y. - 0.25 I X^)/n = -1.35 


5 



0.081 


• Statistics on D (See Appendix B for definitions of terms) 

FC A 


Median, = 10 = 10 

FL 


-1.35 


COV 


= 0.0447 


0^ = g(a,n)s 


g(.01,6) = 1.56 


o = (1.56)(0.081) = 0.126 
o 


Then 




126) /. 434 


<\j 


-1 = 0.30 = 


0.30 




that the R-F method does an adequate job of estimating reliabilities 
at computer costs far leP/S than Monte Carlo, A summary of the R-F 
method is given in Appendix C, 

3) The Wu Algorithm . This scheme was developed by Y.-T. Wu on the same 

NASA/Lewis grant which sponsored this project [10], This numerical 

method, summarized in Appendix D, is more complicated than R-F, 

but increased computer costs, relative to R-F, are insignificant. 

All evidence seems to indicate that the accuracy in estimating 

probabilities of failure is substantially better. At this time 

the Wu algorithm has not been subjected to peer review and has not 

been published, but its performance has been demonstrated to be 

of consistently high quality in a large number of examples. This 

SRP problem is another example. 

The Results . The output of the R-F and the Wu programs are provided 
in Tables 6 and 7 respectively. Results are summarized in Table 5. Agree- 
ment of the three methods in this example is better than usual [10]. In 
this example, three approximations to B and p^ are being compared, although 
Monte Carlo is exact as n ->• «>. 

A practical limitation to Monte Carlo for structural risk problems 
is demonstrated by this example. Note the relatively broad range of 
the 98% confidence interval for a sample of n = 100,000. This range 
would be even broader for the ’’more typical” risk levels of 10 or lower. 

To sharpen the limits, a much larger sample would be required. But even 
for this problem, approximate relative computer costs presented in Table 
8 illustrate the inefficiency of Monte Carlo. 
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In summary, for a single problem, Monte Carlo computer costs may not be 
excessive. But for a large scale program, the R-F and Wu schemes may be 
much more efficient. Furthermore, these methods provide a basis for 
developing safety check expressions for design criteria documents (See 
Example 3 below.) 
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Table 6 


0RIG8MAL tnXJ l3 
OF POOR QUAUW 


Output of Kackwitz-Fiessler Program for Example 1 


COMPUTATION OF THE SAFETY INDEX USING R-F ALGORITHM 


NUMEER OF DESIGN UAR.IABLESi N« 6 
STOP SENSITIVITY = . 00010000 
INITIAL GUESS OF REDUCED VARIABLES 
INITIAL STEP SIZE 0. 1 


0 . 0 


LIMIT 

STATE G(R, S)-0. 

STRAIN RANGE 

PARTITIONING MODEL 



DESIGN 

VARIABLES 


VARIABLE 

TRANSFORMATION 

MEAN/MEDIAN 

COV 

ES 


EVD 

. 53400E-02 

. 20000E+00 

T 


NORMAL 

. lOOOOE+03 

. 50000E-01 

BPP 


LOG 

. 2i600E-0i 

, 99000E-01 

DCP 


LOG 

. 44700E-01 

, 30000E+00 

A1 


LOG 

. 2S100E-'01 

. 47000E+00 

A2 


LOG 

. 15i£»0OE-01 

. &9000E-»'00 

NOTE: 

THE 

MEDIAN IS SPECIFIED 1 

FOR A LOGNORMAL VARIABLE ONLY. 


DESIGN POINT 


VARIABLE REDUCED VALUE BASIC VALUE 


ES 

XR(1)« 

2. 34338 

X(l)« 

. 87259E-02 

T 

XR(2)* 

. 01660 

X(2)« 

. 10008E-I-03 

BPP 

XR(3>- 

78652 

X<3)« 

. 19985E-01 

DCP 

XR<4)« 

. 36855 

X(4>* 

. 49814E-01 

A1 

)CR(5>« 

28086 

X(5>« 

. 24782E-01 

A2 

XR<6>« 

07246 

X<6)« 

. 14907E-01 

SAFETY INDEX . 

BETA 

« 2. 5164 



PROBABILITY OF 

FAILURE 

= . 59346E«02* 




*P^ = 1'(-3) 
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Table 7 

Output of Program Which Uses the Wu Algorithm; Example 1 


COHPUTATION 

OF the SAFETY INDEX 

USING THE Y, T. 

WU ALGORITHM 

LIMIT STATE 0(R 

.S)«0. STRAIN RANGE PARTITIONING MODEL 

DESIGN VARIABLES 

VARIABLE 

TRANSFORMATION 

MEAN/NEDIAN 

COV 

ES 

EVD 

. 53400E-02 

. 20000E+00 

T 

NORMAL 

. lOOOOE+03 

. 50000E-01 

BPP 

LOG 

..21600E-01 

. 99000E-01 

DCP 

LOG 

. 44700E-01 

. 30000E-»-00 

A1 

LOG 

. 2B100E-01 

. 47000E+00 

A2 

LOG 

. 15600E-01 

. 69000E-t-00 

NOTE; THE MEDIAN IS SPECIFIED FOR A LOGNORMAL VARIABLE ONLY. 

DESIGN POINT 

VARIABLE 

REDUCED VALUE 

BASIC VALUE 

ES 

XR(1)» 6.16752 

X(l)*= 

. B6154E-02 

T 

XR(2)» .01855 

X(2)« 

. 10009E+03 

BPP 

XR<3)« -.05075 

X(3)« 

. 19755E-01 

DCP 

XR(4)« .40547 

X<4)« 

. 50354E-01 

Ai 

XR<5)* -.31582 

X(5)-: 

. 24395E-01 

A2 

XR(6)» -.07969 

X(6)*= 

. 14840E-01 


SAFETY INDEX , BETA = 2. 5050 
PROBABILITY OF FAILURE « .&1295E-02* 

*P^ = f(-3) 


A3 


Table 8 


Comparison of Reliability Analyses Between Monte Carlo, 
Rackwitz-Fiessler, and the Wu Algorithms 


Rackwitz-Fiessler 


(a) 


Safety 
Index, B 

2,516 


Probability 
of Failure^^^ 

5.93E-3 


Relative 

Cost 

Factor 

1 


Wu 


(a) 


2.505 


6.13E-3 


Monte Carlo 


(b) 


2.485 


6.48E-3 

98?^ Confidence Limits 
(5.90, 7.10E-3) 
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computed first. Then where $ is the standard 

normal distribution function 

^^^p^ computed by counting the number of failures in a sample 
of n = 100,000. Then B = . 
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EXAMPLE 2, Details of constructing the random variable A are described 
in Appendix A for the two variable case and in Appendix B for a single 
variable, The function g(a,n) is introduced to quantify the statistical 
uncertainty component, essentially by enlarging the sample standard 
deviation, s. It is s which quantifies variability in material proper- 
ties, A summary is provided in Table 9. 

As n becomes larger, this statistical uncertainty becomes smaller, 
and g(a, n) approaches unity. In this example, it is assumed that 
the data of Example 1 is now based on large samples so that all statis- 
tical uncertainty disappears (l.e., g(a, n) = 1) . Table 9 summarizes 
those variables and their COV’s in Example 1 for which this error term 
was included. Also shown is the reduction in COV if the sample size were 
large and g = 1 , . . assuming the same statistics for all variables. 

The goal of this exercise is to demonstrate the effect of statistical 
uncertainty on the design. How important is it to the overall reliability 
analysis to increase sample sizes to reduce this statistical error? 

To accomplish this goal (1) the statistical uncertainty component was 
removed from the COV’s of the four variables considered in Table 9, (2) 
the mean value of service strain, ASg was Increased so that the safety 
index was the same as in Example 1, (for both R-F and Wu) . 

The results using both schemes were identical, as summarized in Table 
10. As the statistical uncertainty is removed, the mean value of ASg can 
be increased at the same level of risk. But the increase in allowable 
strain (mean) is only 2.6%. Thus, in this example at least, it seems that 
while statistical uncertainty may strongly influence the COV of a design 
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Table 9 

A Summary of the Effect of Statistical Uncertainty on the COV's 
of the Random Design Factors of Example 1 

For the random variable, X, the coefficient of variation (COV) Is 



where Is the equivalent standard deviation 
^0 = g s 

s = sample standard deviation 

g = factor to account for statistical 
uncertainty in estimating parameters 


Variable 

COV, including 
statistical 

uncertainty, 1 .e, , g > 1 

COV, assume n large 
enough so that there 
is no statistical 

uncertainty; g = 1 

Bpp 

9 . 9 % 

7.0?^ 

Dpc 

30 , 

19. 


47. 

33, 

-2 

69. 

41. 
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factor, its impact on the overall reliability seems almost negligible. On 
the basis of this one example, it would, of course, be dangerous to con- 
clude that statistical uncertainty is unimportant and that small samples 
are OK, Clearly more studies need to be made on this problem. 



Change in the Mean Value of Sr»"vice Stnain, Acg which would be 
allowed at the same level of risk as Example 1. When the Stati 
tical Uncertainty Component of Dp^, , and Ag is Removed. 
(Essentially same results for both Rackwitz-Fiessler and Wu 
algorithms. ) 


Including Statistical Uncertainty (Example 1) 

5.34E-3 

Excluding Statistical Uncertainty 

5.48E-2 

Percentage Increase: Would be equal to the 
percent increase in the requirement for the 
cross sectional area of a tension member if 
stress and strain were linearly proportional 

2.6 
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example 3f This example demonstrates the use of the R-F and Mu al^^orithms 
to derive a safety check expression which could be used for a design criteria 
document. The problem of Example 1 is used. The general method for 
deriving partial safety factors, (PSF's) is described In the literature 
[8, 10, 12], A simple tutorial in PSF's is provided in Appendix E for readers 
who are unfamiliar with the concepts. 

In conventional design practice, typically a single safety factor is 
employed to account for all uncertainty. A more refined criterion, could 
he developed by applying safety factors to each random design factor. 

In theory, a criterion using these PSF's would produce a more efficient 
design. Described in Appendix E is how probabilistic design methods, 
namely the R-F scheme, can be used to derive the PSF-s. 

It should be noted that a probability based design criterion could 
require that the designer compute p^ (or 8) for the component in question. 

The component would be safe if p^ < p^ (or 8 > 8^) where p^ and 8^ are the 
target risk and safety index respi?.:tively . To require a designer to exercise 
skills in probability theory may be impractical. The much more familiar 
format, a deterministic inequality Involving safety factors, is easy to 
understand and use. In summary, a reliability based safety check expression 
is derived, having a format, familiar to designers, such that probabilistic 
and statistical analyses are invisible. 

In this example, the problem is defined as follows: 

1. The limit state is defined by Eq. 17 (with the substitution of 


Eqs. 18 and 19). 
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2. Distributional forms and statistics of the design factors are 
defined in Table 2. 

3. The target safety index, is given as 3^ = 3.0, Just for ref- 
erence, , , e the notional probability of failure associated with this 
value is “ <I> (-3) « 0.0013. 

4. The tuning factor A (See Appendix E) is defined by replacing 
Ae with A»Ae . In the R-F algorithm, A is adjusted so that 3 = 3.0. 

b D 

As an alternative viewpoint, Ae can be taken as unity, and A then could 

D 

be thought of as the mean of Ae . 

5. The nominal values are defined as the median values for variables 
having lognormal distributions and the mean values for the other variables. 

With this information, the PSF program at the University of Arizona 
was run; the results are presented in Table 11. Input to the program are 

the variables, their distributions and statistics, the nominal values 

(at the top of the table) , and 3^ (at the bottom) . The program computes 
the partial safety factors, as listed. 

Combining the PSF’s with the limit state expression, a safe design 
results when the following inequality is satisfied. 

1.84 (M ) < {0.91 B exp l-1.14Dp. f*]) ( 21 ) 

1.0. 86A, O. 95 A 2 J 

where , 

fpp =1.1-0.10 log^Q T 

f = 1 - f 
PC PP 

Note that the relationship for f^p is valid only for a limited range of 
hold time (See Fig. 12). 
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ORIGINAL PACaS iiS 
OF POOR QUALITY 


Table 11 

Results of PSF Program 


FAILURE FUNCTTQ^J: SRP MODEL 


OESIGM VARIABLES 


VAR lABLF 

OTSTRIBUTIDN NOMINAL 

MEAN/MEDIAN 

COV 

ES 

EVD 

.lOOOOE+01 

.IOOOOE4-01 

.20000 

T 

NORMAL 

.10000E*03 

.lOOOOE+03 

.05000 

BPP 

LOG 

.21600E-01 

.2160CE“0;i 

.09900 

OCP 

LOG 

.^4700E-01 

.^<r700E-01 

.30000 

A1 

LOG 

.28100E-01 

.28100E-01 

.A7000 

A2 

LOG 

. 15t)00E-01 

.15600E-01 

.69000 

NOTF: the 

MEDIAN IS 

SPECI FIED'^TDP A 

LOGNORMAL VARIABLE 

ONLY. 


(y) 


VARIABIF 

REO'JCEO 

3FSI GN 
VALUE 

POINT 

BASIC VALUE 

PARTIAL 

SAFETY FACTOR* 

F. S 

XR(l) . 

2.B0029 

X(D- 

.18404E+01 

1.8404 

T 

XR(?) - 

. 02031 

X( 2 )* 

.lOOlOE+03 

1.0010 

BPP 

XR(3)« 

-.92031 

X( 3)- 

.19723E-01 

.9131 

DCP 

XR(<»I ■ 

.A^062 

X(4I- 

.50B73E-01 

1.1381 

A1 

XR (51 - 

33125 

X (5)* 

.24235E-01 

.6624 

A2 

XR(6>- 

-.08437 

X { 6 ) * 

-14785E-01 

.9470 


safety index# beta* 3.00 
SCA'LF FACTOR- .A6658E-02 A 


♦Assumes that nominal values = 


0 

median for lognormal variates 
mean for Ae^ and T 
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where the bar over a variable denotes mean, and the tilde denotes median 
(the nominal values by definition) . 

This example was presented only for demonstration purposes . Limitations 
of its use in a design criteria document center around the fact that each 
PSF is a function of all of the statistics and parameters. 

1. The expression was derived on the basis of known statistics of 
the design factors. It applies to a specific case* For aoplication, a 
range of possible statistics and corresponding PSF's should be studied 

to construct characteristic PSF's. The problem may require some engineering 
judgement in smoothing the PSF's. 

2. The PSF's were derived for a specific life, N = 1,000 cycles. 

If the requirement should include other values of N, then the behavior 
of the PSF's should again be scrutinized. 
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SUMMARY COMMENTS 

Reliability technology has now developed to the point where application 
to complicated problems is a practical reality in many cases. The Rackwitz- 
Fiessler and Wu algorithms provide an estimate of reliability of a component 
experiencing the combined effects of creep and fatigue. The strain range 
partitioning form of the limit state has a relatively complex and highly 
non-linear form; yet, as demonstrated herein, these algorithms easily handle 
this problem with negligible computer costs. 

Reliability analysis was used to assess the impact of small sample sizes 
on component risk. In addition to uncertainty due to inherent data scatter 
resulting from material behavior, statistical uncertainty, resulting from the 
fact that parameter estimates are random variables, is present. An example 
provides an illustration that statistical uncertainty may be relatively 
insignificant, but it would be dangerous at this time to present this as 
a general conclusion. 

These advanced reliability methods can also be used to derive safety 
check expressions which employ partial safety factors. A maximum allowable 
risk is the basic criterion. An example in which PSF's are derived for the 
SRP problem was presented. As a general comment, caution should be exercised 
in specifying PSF's for general application simply because they are functions 
of all of the statistical parameters in a given limit state. 
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APPENDIX A. ANALYSIS OF STRAIN-LIFE DATA 

Methods for analyzing strain-life fatigue data are discussed in this 
appendix. The goal of such data analysis is to provide a characterization 
or statistical suininary of the e-N relationship in a form suitable for inclu- 
sion in a comprehensive reliability analysis. 

Considered will be uncertainty associated with (a) inherent behavior 
of the material as evidenced by scatter in the data, and (b) statistical 
behavior of the least squares estimators. This problem was addressed in 
Ref. 11, and the following summary describes a model for quantifying both 
uncertainties, thereby producing a model for reliability analysis. 

The Least Squares Line 

Consider a constant amplitude fatigue test in which pairs of data 
(Ae^, N^) , i = 1, n are collected. N^, is the cycles to failure associated 
with strain-range (or amplitude , Ae^) and n is the sample size. Ae is 
the independent (or controlled) variable and N is the dependent variable. 
Hypothetical test data are shown in Fig. A.l plotted on log-log paper. 

There data imply a model of the form 

N = A(Ae)^ (A.l) 

where A is a random variable and b is constant. Therefore, N would be 
a random variable also; its density function is shown in Fig, A.l. 

Because b and A would be those parameters in a design algorithm which 
represents the fatigue strength of a material, it is necessary to provide 
a description of b and A. 
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Therefore, the problem is to translate the data , N^) into the 
value of b and a distribution of A, In order to do this, first consider 
the median curve defined by A (the median of A) in Fig. A.l. A linear form 
of this median curve is 

Y = a + bX (A. 2) 


where 

Y « log N, X = log Ae, a = log A (A. 3) 

Eq. A. 3 translates the data N^) into (X^, Y^), i = 1, n, Eq. A, 2 

defines the mean of Y (log N) given X (log Ae) . The scatter in the data 
is defined by the standard deviation of Y given X, denoted as a and assumed 
to be constant (not a function of X). Moreover, Y is assumed to have a 
normal distribution for all X. 

A 

Using the method of least squares, a, b, and a are estimated by A, b, 
and § respectively [11], 

b = ? (X - X)(Y - Y)/ I (X - X)2 
i=l ^ i=l ^ 


a = Y - bX 


2 

s 




Ca + bX.)]^ 


(A. A) 


where X and Y are the sample means of X and Y respectively. Because each 
Y^ is a random variable, the estimates d, b, and s are also random variables. 
The "best fit” line 


Y = a + bX 


(A. 5) 
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is called the least squares line, Y is the estimate of the mean of Y given 

X. 

Note that (a) Y given X is normal and (b) the least squares line is 
the estimate of the mean of y|x. Therefore, it follows that (a) NjAe 
is lognormal and (b) the least squares line, N, is the estimate of the 

A 

median of NjAe. As illut 'rated in Fig, A. 2, the least squares line, Y, 

A 

is only an estimate of the actual median by virtue of the fact that a 

A 

and b are only estimates. 

The general goal of this study is to develop an empirical relation- 
ship between Y and X which accounts for both the scatter in the data and 
the distribution of the estimators, but is easy to use in probability- 
based design formats, A proposed model, suggested by the above discussion, 
is as follows: 

A 

1. Let b = b be a constant. 

2. Assume that the uncertainty due to both sources is accounted for in 
a (and therefore A) by considering the y intercept as a random 
variable, 

3. Therefore, let the empirical relationship be 

Y == a + bx (A. 6) 

o 

where a has a normal distribution with mean a and standard 
o 

deviation a . 

o 

The concept of an equivalent prediction interval (EPI) was employed 
to derive a^[ll]. 
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“ 8(n,a)s (A. 7) 

where 

g(n»Oi) » exp[A(a){ln 

A(a) - l,56[i (A. 8) 

B(a) - 3.32 - 1.7a 

6 ;ln<50; 0.01 iSci £0.15 

g(n,a) is in essence, an adjustment factor to s to account for the fact 
that there is uncertainty in the pstlmates of a and b and s. In turn, 
s accounts for the scatter in material behavior. 

The value of a is arbitrary. It refers to that region of the tall 
area where it is desired to have a good fit [ 5, ’1]. As a general rule, 

a value of a » 0,01 is reasonable. For reference, g(a,n) is plotted as a 

function of n for a ® 0.01 in Fig. A. 3. g is the measure of statistical 

error, and it is interesting to note how quickly it drops as n increases, . 

thus suggesting that statistical uncertainty may be small for n > 10, 

The consequences of the model described above, relative to reliability 
analyses are: 

1. y|x has a normal distribution. (Thus N given As has a lognormal 
distribution) 

A 

2. The mean value of Y|x is a + bx. (Thus the median of N is 

A 

= 10® (Ae)^ 


3 . 


The standard deviation of Y X is a (and is not a function of X) . 

o 
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4i " log A is normal and A is lognormal. The median % and 
coefficient of variation of A can be obtained from the log- 


*H>rraal (base 10) forms 


X » lOa 


(A.9) 
(A. 10) 


EXAMPLE Given the PP strain-life data (n * 9) as illustrated in Fig. 13 
and given in Table 3 , it is required to produce to statistics on A. 

From the least squares analysis , 

a “ -1.552 b - -1.57 s - 0.138 

The median of A is computed by Eq. A.9, 

f\j -1 552 
^ = 10 « 0.0281 

From Fig. A. 3, for a = .01 and n « 9 

g(.01, 9) = 1.41 

The equivalent standard deviation is, 

Go = 8 • S 

= (1.41) (0.138) = 0.194 
and the COV of A is computed from Eq, A. 10 as, 

c = 4 o0.194^/.A34 = 0.472 
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APPENDIX B* STATISTICAL ANALYSIS OR SINGLE RANDOM VARIABLE 

ConSju,i«!r a random sample of size n of a single random variable, X. 

I « (X^, .1^, • * • V 

It is known that X has a normal distribution. The sample mean, X, and 
sample standard deviation, » , are computed, To establish a "design 

X • 

value" of X, the notion of an equivalent prediction interval (EPI) can 
be used [ 5 ] , But the EPI can also be used to provide "improved" statistics 
for probabilistic design, 

Define an equivalent standard deviat.** n as, 


where 




'a;n-l 


A + o7n) 


t , = students t variate 
a;n-l 


z = standard normal variate 
a 


(B.3) 


a = reference probability level 

The choice of a is arbitrary, but for general design, a value of 0.01 is 
recommended. Reference 5 provides additional discussion. For conveni- 


ent reference, the value of gj|^for a = 0.01 is presented as a function of 
n in Fig. B.l. 
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For design purposes, one can state that X has a normal distribution 
having mean and standard deviation (X> a^) . Scattf :.i "''erent in the 
phenomena is described by s^, and statistical scatter, i.e., the fact that 
X and a are random varlableo is described by g. 

Example: Let Z be a lognormal ly distributed random variable. Let 

Y ® log^QZ, Then Y has a normal distribution, A random sample, of 

size n ® 6 is taken. Transforming to Y and computing the statistics, 

Y = -1.807 and s = 0.173, 

The equivalent standard deviation of Y is given by Eq. B.2, From 
Fig. B.l, gj^= 1.56 for a «= ,01 and n = 6. Thus, 


a = (1.56)(0.173) 
o 

= 0.270 

•Invoking basic properties of the lognormal distribution (See Appendix C) . 


• Median of Z 

^ = 10^ = 0.0156 

• Coefficient of variation of Z 


« / 10 ^ 0 ^* 


434 


« 0.687 


-1 
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APPENDIX C. THE RACKWITZ-FIESSLER (R-F) ALGORITHM 

The algorithm proposed by Rackwitz and Fiessler (7) has been ejctensively 
described in recent literature (8, 9, 12). The procedure for calculating 
the R-F algorithm safety index can be summarized as follows: 

1. Define each design factor, X^^ (i»l,n) and its corresponding 
probability distribution; F^ and f^ denotes to cdf and pdf of 

respectively. 

2. Define reduced variables 


u, = i i 
i 


i = l,n 


(C.l) 


where = mean and standard deviation of X^ respectively. 

3. Define the limit state in reduced variables 


g’(u) = 0 


(C.2) 


where u = (u , u_, . . . u ) 
'VIZ n 


4. Make an initial estimate of the safety index 


• 

an vu. 


rain /u^ 


. + u‘ 


(C.3) 


subject to This is the Hasofer-Lind generalized safety index. 

5. Calculate the corresponding design point, X where 


X , = u . 0 , + p . 
X i i i 


X ** 1 y III 


(C.4) 


The design point is that point on the failure surface closest to 
the origin of y educed coordinates. 
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ORIGINAL PmE ^ 
OF POOR QUALITY 


6. Calculate the means and standard deviations of the equivalent 
normal distributions for each non-normal variable 


a 


Ni 



(Fj,<X*))] 




(C.5) 


where <j) ® standard normal pdf and ^ - standard normal cdf 
for each variable. 

7. Define the new reduced variables 


I 




(C.6) 


8. Calculate a new estimate of the safety index 

3^ = min/(u|)^ + + . . . + (u^)^ (C.7) 

subject to g’(u’) - 0. 

9. Repeat steps 6 through 10 until the difference 

Ib - S I < t (C.8) 

where t is the "error." In this study, the value t = 0.001 
was used. 

A 

10. The probability of failure Is calculated using 3 = 


<I>(-3) 


(C.9) 
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APPENDIX D. A NEW METHOD OF CONSTRUCTING EQUIVALENT NORMAL DISTRIBUTION 
AND COMPUTING PROBABILITY OF FAILURE 

This method constructs, for each non-normal variable, a scaled 
three parameter equivalent normal distribution function (cdf) employing 
a least square scheme. The probability of failure p^ is computed using 
the obtained parameters by assuming that the limit state is linear at 
the design point. 

Consider a limit state function g(X) involving n Independent random 
variables which is linearized at the design point, 

gC5) = + (D.l) 

Three equivalent normal parameters, (A, Ojq) » each X^ are found, 
one by one, and the probability of failure is estimated as 

Pf " f " f ‘^u ^i^ 

^ 1=1 

where A. are the scale factors, fy(;i<) and (J) (u) denote the joint prob- 
^ ^ ''si 

ability density function (pdf) of the original and equivalent normal 

region on the original 

failure region on the 


variables, respectively, is the failure 

coordinates, and ^^2 corresponding 

reduced coordinates, u, in which 

’ r\j 


Pi 


^1 ~ ^Ni 
°Ni 


(D,3) 
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() 


Using Eq. D.3, the reduced limit state function is also linear. A minimum 
distance, 5, on the u space can be found, and is estimated as 

n 

p B $(-6) (IIA ) (D.4) 

1-1 

The generalized safety index, 3, is computed as 




(D.5) 


Consider constructing the equivalent normal cdf for one of the 
variables, X^, and let 


Y = a.X, + . . . + a . - X. , + a.,, X, + , , . a X 
1 1 i-1 i-1 i+1 i+1 n n 


(D.6) 


Eq, D.l becomes 

g(J) = + Y 


(D.7) 


Thus, the limit state involves only two variables; X^ is the variable to 
be normalized and Y represents the sum of the other variables. Assume 
that the Rackwitz-Flessler (R-F) algorithm has been performed, and ^ are 
replaced by the equivalent normals , then Y is also a normal variable with 
pdf of 





(D.8) 


which will be used in the following procedure. 
Define the R-F reduced design point as 

A 



(D.9) 
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where and are the R-F equivalent normal mean and standard deviation, 

respectively. A variable may be defined as a ’’strength" variable if < 0, 

and a "stress" variable if z. > 0. 

1 

Assume that a non-normal variable, denoted as X (without the sub- 
script i) is a strength variable with cdf of F(x) and pdf of f(x). The 
three equivalent normal parameters can be found by minimizing the sum 
of the errors of the squares between two functions, F(x)(|i’(y) and 
A$(x)(})’(y), l.e., 

Min: E « / [A<Kx)(J) ’ (y) - F(x)<J> '.(y) ]^dx (D.IO) 

Subject to g(^) «a^+aX+Y«0 


where A<I>(x) is the equivalent normal cdf with mean and standard 
deviation 

The procedure described in the follolwng imposes two constraints, 
similar to the R-F algorithm, to Eq. D. 10, i.e., match cdf’s and pdf’s 
at the design point , 

A<ji(u*"^) «= F(x ) (D.ll) 




f(x*) 


(D.12) 


where 4»(u ) and <{>(u ) are the "standardized" normal cdf and pdf, respectively. 
Using Eq. D.ll and Eq. D.12, the error sum, E, can be evaluated for a 
given A value. 
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The procedure 
as follows: 

1. Calculate 
Define 


for determining the normal parameters can be summarized 


as a function of x. 


z 


1 



(D.13) 


z 


2 



(D.IA) 


The reduced limit state using R-F results can be derived as 


so that 


* * 2 
^1 ^ 



(D.15) 


(D.16) 


Given any x value, z^ can be calculated. Therefore, <}>'(y) 
can be computed using Eq. D,8, Note that because o^, in Eq. D.8, 
is a constant, it can be taken out from the E integral without 
affecting the result of the parameters. 

2. Make an initial guess of A (e.g., A = 1) and calculate u from 


Eq. D.ll, 


* 


u 




3. Calculate a„ from Eq. D.12 
N 


a 


N 


A 

f(x ) 


(D.17) 


(D.18) 
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4. Calculate using Eq. D.3, 

]x^ • Y. - u0j^ (D.19) 

5. Compute the error sum, E, in which 4>(x) is evaluated using 

and a... 

N 

6. Choose other values of A and repeat step 2 through step 5. An 
"optimum" A, which minimizes E, can be determined using a suitable 
optimization routine, Three parameters are thereby determined. 

7. Repeat the above procedure for other non-normal variables. If the 
variable is a stress variable, E(x) and $(x) should be replaced by 
1 - F(x) and 1 - 'I'(x), respectively, in all the formulations. 

8. Compute p^ and 3 according to Eq. D.4 and Eq. D.5, 

Because, in general, there is no closed form solution for the E 
Integral, a numerical scheme must be used to approximate E by replacing 
the integral by a summation and replacing dx by Ax, i.e., x values must 
be discretized. The region of x must also be set. It can be determined 
such that 

< H for the two limits of x (0,20) 

F(x ) 

where H is a reasonably small value, say, 0.2. Note that F(x)<j)’(y) relates 
closely to p^, therefore small H value implies that a sufficiently wide 
region of x will be Included in the summation of E. However, when a variable 
has a relatively small z (e.g., z^/3 < 0.1) and a large coefficient of vari- 
ation (e.g., 0.4), the range of x may become very wide (therefore, too many 
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points of X need to be included in the summation) to satisfy Eq, D.20. 

In such cases, it is suggested that R-F equivalent normal pareuneters may 

A 

be used (i.e,, A ■ 1) directly. Because is small, the difference in 
p^ estimate is usually negligible, 

A user-oriented computer program applying the above numerical scheme 
has been developed in the University of Arizona, The process of choosing 
X values is automated; only the distributional information and the limit 
state need be input by the user to generate the probability of failure 
estimate. 
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APPENDIX E. A SIMPLE ILLUSTRATION OF THE COMPUTATION OF PARTIAL SAFETY 
FACTORS 

Definition of and derivation of partial safety factors (PSF) is described 
in the literature [8 » 10»12]* Because techniques of constructing design 
criteria using PSF's are relatively new and not widely known, in the aerospace 
and mechanical design community, a simple example is presented which provides 
a tutorial. 

THE EXAMPLE 

The limit state is 


L + D « AR 


where L and D are stress variables and AR is strength variable. 
The target safety index is chosen as 

3 “3.0 
0 


(E.l) 


(E.2) 


The statistical parameters are given as follows; (y “ mean, o “ standard 
deviation, tilde indicates median, C = COV) 


L 

Extrene Value 
Distribution (EVD) 

'^L * 

'^L " 

2.0 

D 

Lognormal (LN) 

5 “ 20 

S ■ 

0.15 

R 

Weibull (WEI) or 
Lognormal (LN) 

Mr - 50 

i 

°R “ 

5.0 


A = constant; here we could assume that it is s geometric variable, 
e.g, , cross sectional area. But in the process of computing 
the point, A plays the role of an adjustment or "tuning factor." 
Its role will be described later. 



73 


now THE PARTIAL SAFETY FACTORS ARE COMPUTED 

The Rackwit 2 -Fics 8 ler (R-F) method for approximating non-normal variates 
with an equivalent normal in a Hasofer-Lind generalized safety index approach 
is employed. An R-F program, (named RACA) , based on an optimation method 
for computing the safety index, S, has been developed at the University of 
Arizona. This program was used for the calculations, ^SF^s v?ere com- 

puted by the following steps. 

First RACA was used to compute the design point so that $ ■ 3.0, the 
target safety index. This has to be done (with the present version of the 
program) by iteration by adjusting the value of A so that 3 * 3.0. Thus, A 
is called the tuning factor . The output of the program is Illustrated in 
Table 2. The results of the program are the design point (L>'f, D*, R*), , . . 
and A. Note that the design point is on the failure surface, i.e., 

"stress" "strength" 

Define partial safety factors , 

L* » D* = (AR*) » 

where the subscript "n" refers to the nominal value. This value is ar- 
bitrarily chosen. It could be chosen as the mean or median, or perhaps 
a value in right tail for stress variable or in left tail for strength 
variable. Clearly the partial safety factors depend upon the definition 
of nominal values and therefore, codified safety check expressions should 
clearly specify the definition of a nominal value. 
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Substituting Eq, E,A into E,3 
But if m Xet 

Y-E + Yn^ ^ Ytj(AR ) (E,6) 

VL n 'D n - 'R n 

we will Insure that 0 > Here we are just saying that it is okay to 
lower stresses or increase strength. 

In this example, A » 1.06, and 



L* 

* 11.17 


Design Point 

D* 

•“ 23.53 

(E.7) 


AR* 

= (1.06)(32.73) 


Assume that means 

are nominal 

values (often the mean val,. 

Is used 


for stress variables, but some number in lower tall for strength variable). 
Then it follows from Eq, E.4 that 


L* 

11.17 - 


10.0 

D* 

23.53 . 


20.0 

AR* 

. 32.73 


50.00 


1.12 


1.18 


- 0.655 


Thus, the safety check expression or condition for a eafe design Eq. E.5 
expressed in terms of the PSF's, becomes, 
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1.12 L + 1.18 D < 0/65- (E.9) 

n n “ n 

Eq, E.9 could then be employed as a. safety check expression in a 
design criteria document, It is of course understood chat the inequality 
would be valid only for cases when the design factors are assumed to possess 
the same statistics as the variables used to derive the PSF’s, 
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